I will first tidy the data set. It is important to standardize some of the values including, snow, tmax, tmin, and prcp.
snow
tmax
tmin
prcp
ny_noaa %>% plot_ly(x=~year, y = ~snow_mm, color = ~year, type = "box", colors = "viridis")